Non-Linear a Model for Inhomogeneous Spin Chains 
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We derive a non-linear a model (NLSM) generally representing antiferromagnetic Heisenberg spin chains 
with inhomogeneous spin magnitudes and inhomogeneous nearest-neighbor exchange constants arrayed in 
finite periods. Only a restriction is that the average of spin magnitudes on one sublattice is the same as that 
on the other sublattice. The NLSM gives the gapless condition explicitly written by the spin magnitudes 
and the exchange constants. We apply this gapless condition to several cases including systems with 
impurities. 

PACS numbers: 75.10.Jm, 75.30. Et, 75.30. Hx 



In 1983 Haldane predicted that an antiferromagnetic 
Heisenberg spin chain is gapful if the spin magnitude is an 
integer, and is gapless if it is a half-odd-integer [Q. Since 
his theory is based on mapping of a spin model to the 
non-linear a model (NLSM) , the usefulness of the NLSM 
is recognized by researchers in condensed matter physics. 
Various interesting aspects of the method of the NLSM 
are mentioned in some books Extension of the 

NLSM method to systems with inhomogeneity in spins 
and exchange constants is a challenging problem. It is not 
only theoretically important but also useful in explain- 
ing and suggesting experiments. Affleck reformulated the 
NLSM method and applied to a spin chain with bond al- 
ternation |§. He divided the spin chain into blocks each 
of which contains two spins, and defined new variables for 
each block. Affleck and Haldane analyzed and predicted 
gapless points for general homogeneous spin chains ||. 
Fukui and Kawakami tried to extend the NLSM method 
to a spin chain including impurities and to spin chains 
with unit cell containing three and four spins . They 
rather intuitively introduced new variables to form an 
NLSM and did not care to preserve the degrees of free- 
dom. Hence we cannot determine whether their resultant 
NLSM's are correct or not by inspecting their theory it- 
self. 

In this Letter, we obtain a NLSM generally represent- 
ing a wide class of inhomogeneous antiferromagnetic spin 
chains. The derivation is based on the idea dividing the 
system into blocks || but blocks used here are generally 
defined to contain more than two spins. In a path inte- 
gral formula, the variables of integration for each block 
are transformed to form the NLSM with preserving the 
original degrees of freedom. We apply this NLSM to sev- 
eral cases including systems with impurities and examine 
gapless conditions. 

We consider the spin chain represented by the Hamil- 
tonian 



N 



(i) 



where Sj is the spin at site j and Jj(> 0) is the exchange 
constant between Sj and Sj+i. The number of lattice 
sites is N, the lattice spacing is a and the system size is 
L = aN. The quantum number for the magnitude of Sj 
is denoted by Sj. The system is periodic with period 2b 
(b: a positive integer): 



Jj+2b — Jj, Sj +2 b — Sj. 



(2) 



We divide the spin chain into blocks each of which con- 
tains 2b spin sites. The block size 26a is a positive integer 
(even number) times the size of the unit cell if a unit cell 
contains an even (odd) number of spins. In a block {Jj} 
and {sj} are arbitrary except for the restriction: 



S 2 j-1- 



(3) 



This restriction excludes systems with ferrimagnetic 
ground states in the classical limit. 

The partition function Z is written in a path integral 
formula. When (— l) 3 Sjiij is the expectation value of Sj 
for a coherent state, it is 



Z = J D[n 3 ] 115(^-1) 



(4) 



The action S at temperature 1/(3 is given by 

/n 
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S = iJ^-lJV^] + - / oV]T J, (<5n,) 2 (5) 

3=1 J ° 3=1 



3=1 

with Srij = nj + i — n-,- and 



1 



Jj = JjSjS j+ i- (6) 

The first term in Eq. (||) comes from the Berry phase and 
w[rij] is the solid angle which the unit vector rij forms in 
the period [3. 

We transform spin variables {n, } into gradually chang- 
ing variables {m(p)} and small fluctuations {L q (p)} : 
where p labels a block in the lattice (p = 1, 2, • • • , N/2b) 
and q labels a site in the block (q = 1, 2, • • • , 26). The 
spin variable at the qth site in the pth block is written as 

n 2b P + 9 = (l -7?^^)m(p) +7 9 ^^m(p-7,j) 

+(-l)«oL,(p), " (7) 

where 7, = +1(-1) torl<q<b(b + l<q< 2b). The 
original constraints {n| = 1} are changed to {m 2 (p) = 1} 
and {m(p) • L g (p) = 0}. Here we notice that the num- 
ber of variables increases in this transformation. To solve 
the problem we can add an additional constraint for each 
block which is arbitrary as long as it consists with the 
other constraints; e. g. l) 9 k g (p) = or L 2 b(p) = 

is a possible one. Hence 2b vector variables {*i2bp+q \ q — 
1, 2, • • • , 26} per block are transformed to 2b vector vari- 
ables {m(p), L g (p)| q = 1,2, • • - ,26- 1} per block; L 2 &(p) 
is written by the other variables through the additional 
constraint. Also 26 original constraints {n q = l\q = 
1,2, •••,26} per block are replaced by 26 constraints 
{m 2 (p) = 1, m(p) • L,(p) = 0| q = 1, 2, • • • ,26 - 1} per 
block. Thus we obtained a new set of variables without 
changing the original degrees of freedom in the path in- 
tegral formula (M). In other words, we only transformed 
the variables of integration. As for the additional con- 
straint, the choice does not affect physical quantities at 
all. This is because the variables {L g (f>)| q = 1, 2, • • • , 26} 
only appear as {L g (p) + L q+ i(p)\ q = 1, 2, • • • , 26 — 1} in 
the partition function, as will be seen. 

To take a continuum limit wc identify the center of the 
pth block, (2bp+b)a, as coordinate x. Then the difference 
between adjacent spin variables is replaced as 

Sn 2bp+q -> a[d x m(x) - R g (x)] (8) 

with 

R q (x) = (-mL q (x)+L q+1 (x)}. (9) 

Equation (Q) for q = 26 reads as T{.2b(x) = L 2 &(ir)+Li(a;). 
In the Berry phase term of the action (||), the following 
relation stands due to the restriction (||): 

2b 2b 

^2(-l) q s q w[n2b P + q ] = y]s q dw[n 2bp+q ], (10) 

q=l q=l 

where 5w[n 2 b p +q] = w[n2b p + q +i] - w[n 2 bp+ q ] and 

§ 9 =EH) W ^. (11) 
fc=i 



The continuum limit of the Berry phase term is taken af- 
ter 5w[a.2bp+ q ] is transformed to the r-integral of dn2b P + q - 
(n-2bp+q x d T ri2bp+ q ) as in the usual way §{§. Thus the 
action (|J) becomes S c = Si + S2 with 




where s' = Y%=i Sq/b, J = Y^Li Jq/^b and i q = 
iaJ q (d x m) + s q (m x d T m). 

The variables {L 9 } appear only as {R g } in the action 
S c - The variables {R 9 } are not independent and equa- 
tion Y^ q =i Rg = stands due to the definition (||); e. g. 
we can delete R2b by the equation. We treat the equa- 
tion as a new constraint. The constraints {m • L g = 0} 
are rewritten as {m • R q — 0}. Instead of deleting some 
variables by constraints, we insert the corresponding 8- 
functions into the path integral formula (||) and treats all 
the variables independently. We use the integral repre- 
sentations of #(X)g=i R</) an d ^( m ' Rg) with integration 
variables u and a q . Then the following term appears in 
addition to the action S c : 




Carrying out integrations in the partition function first 
with respect to {R q } and then to u, S2 + S3 reduces to 




where F 9 = f q + a q m and 

2b 2b 

F^t^^/^fi,)- 1 . (16) 

q=l 9 =1 

Expanding the integrand in Eq. (|l"5|), we find that {a q } 
appear only in a bilinear form and are integrated out. 

Collecting Si and the remnant of S' 2 after the inte- 
gration with respect to {a g }, we have the final effective 
action: 

rP rL f j(0) 

S c s = J At J dxi — i-j^-m ■ (d T m x d x m) 

+ ^(^^W + i J(0 '<M' (17) 

where are defined as 
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1 y (~s q ) n 
26 ^ 

q=l 



Ja 



(n = Q,l,2) 



(18) 



with Eqs. (|6j) and (|1 1|) . Thus we have obtained the ac- 
tion of the NLSM describing the Hamiltonian ([!]) in the 
continuum limit. The real-space cutoff is the length of 
a block of the minimum size. The topological angle 9 is 
given by setting j(°)/JW in the first term as O/Att. The 
velocity v and the coupling constant g are given by equat- 
ing the coefficients of (<9 T m) 2 and of (d x m) 2 as l/2gv and 
v/2g respectively. Note that the coefficient of (<9 T m) 2 is 
always positive. 

The action ( fL7| ) is independent of the way to divide the 
system into blocks. First we displace each block by one 
site. Then the spin magnitudes in a new block are or- 
dered as (s 2 , s 3 , ■ ■ • , s 2 b, si) instead of (si, s 2 , s 3 , ■ ■ • , s 2 &). 
Denoting quantities related to the new blocks by let- 



ters with prime, we have relations §' = s 



_i for 



1 < q < 26 — 1 and s' 2b — 0. Using these relations we 
obtain the following transformation: J'W = , 



1 



Ji I L = A _ + A_ n 9) 

J'U) J(0) J(l) ' J'(2) j(0) j(l) + j(2) ■ ^ 1J ' 

This transformation do not change the coefficient of 
(<9 r m) 2 in Eq. (|l7|). The topological angle (divided by 
2-k) changes as 



2J'(°> 
J'(i) 



2s x 



2JW 



(20) 



Here the first term of 2s i is an integer and does not 
affect physics. The negative sign of the second term is 
also irrelevant. Thus the action ( |l7j ) is invariant under 
the block displacement. 

Second we inspect the effect for the action when we 
use the block of size 2rba (r: a positive integer) instead 
of 2ba. The order of the spin magnitudes in a new block 
is then r times repetition of (si, s 2 , ■ • • , s 2 &). Because 
of the restriction (^|) we have relation s' 2 j b+q = s q for 
q = 1, 2, • • • , 26 and j = 1, 2, • • • , r. Using this relation 
we see that J'w = Jw for n — 0, 1,2 and the action 
( |l7| ) is invariant under the block enlargement. 

The NLSM has a gapless excitation when 0/2ir is a 
half-odd-integer. This condition is written as 



2 J(°) _ 21 - 1 



(21) 



where I is an arbitrary integer. In what follows we exam- 
ine this condition for several cases. 

We first apply the general formula (|l7]) to the ho- 
mogeneous case: si = s 2 = ■ ■ ■ = syy = s and 
J\ = J 2 = • • • = Jjv = J- I n this case Eq- Jig ) gives 
j(°) = J S 2 , j(!) = 2Js and J< 2 ) = 2 J. The coefficients 
in Eq. ( |T7j ) are simple as 0/47T = s/2, l/2gv = 1/SaJ and 
u/2g = aJs 2 /2. Hence Eq. ( O ) in this case is equivalent 
to the NLSM which Haldane originally considered 



For 6 = 1, a block contains only two spins and the 
restriction (||) reads as si = s 2 . Although the cor- 
rect NLSM in this case has been already obtained 
we restate some results based on the general formula 
(0). { J(")} are calculated as = 2s 2 J 1 J 2 /(J 1 + J 2 ), 
= 2s Ji and — 2J\. The gapless condition 

@ is then sJ 2 /(Ji + J 2 ) = (22 - l)/4. This equa- 
tion gives the single gapless point J 2 j J\ = 1(1 = 1) for 
s = 1/2 and J 2 /Ji = 1/3(2 = 1) for s = 1. Numeri- 
cal calculations for s = 1 show that the gapless point is 
at J 2 /J\ = 0.6 d,[l(J and an experiment for [{Ni(333- 
tet)(/Li-N 3 )}„](C104)„ agrees with this value Q. Hence 
the method of the NLSM does not always give quantita- 
tively correct results. However the NLSM is expected to 
represent the essence of quantum spin systems. 

For 6 = 2, the restriction (||) is now s\ + S3 = s 2 + S4. 
In the case of s% = s 2 , S3 = S4, J\ = J3 and J 2 = J4, we 
have 2J(°)/J W = 2sis 3 (si+s 3 )/(s 2 + s 2 + 2sis 3 Ji/J 2 ). 
Following the condition (|2l|), a gapless excitation appears 
at J 2 /J\ = 4/7 for s\ = 1/2 and s 3 = 1, and at J 2 /J\ — 
3/7 for si = 1/2 and s 3 = 3/2. For s x = 1 and s 3 = 2, 
two gapless points are J 2 /J\ = 4/19 and 4/3. Fukui 
and Kawakami || have obtained the same results Jl2| . 
Tonegawa et al. performed numerical calculation for s\ — 
1/2 and S3 = 1, and obtained J 2 /Ji — 0.77±0.01 for the 
gapless point [^3|. In the case of si = s 2 , S3 = S4, J\ = J 2 
and J3 = J4 , we have 2 jW/J^) = 2sis 3 /(si+s 3 ). Hence 
the gapless excitation appears if and only if si = S3 and 
the value is a half-odd-integer irrespective of the values 
of J\ and J3. In the case of si = s 2 = S3 = S4 = s, 
J\ = 1 — <5, J3 = 1 + <5 and J 2 = J4 = J, we have 
2J(°)/JW = 2sJ/(J+l-5 2 ). Chen and Hida performed 
numerical calculation for s = 1/2 and obtained a phase 
boundary fll4[ . The positive- J part of their boundary is 
close to S = (1- J) 1/2 determined by Eq. (|l]) with 2 = 1. 

Using the block of 6 = 3, we can deal with various 
systems. Here we examine the case that a unit cell 
contains three sites. In a block of the minimum size, 
the spin magnitudes and the exchange constants are or- 
dered as (si, s 2 , s 3 , si, s 2 , s 3 ) and (Ji, J 2 , J 3 , Ji, J 2 , J3). 
The restriction (0) is satisfied. In this case we have 
2j(°)/j( 1 ) = si — s 2 + S3 irrespective of the values of 
Ji, J 2 and J3. The condition ( ^l|) says that systems 
with one or three half-odd-integer spins in a unit cell are 
gapless 



15 



We examine a bond impurity system in the general for- 
mula (|l7| ) . In this system all the spin magnitudes are the 
same and denoted by s; {s q } are the same as those for 
the homogeneous case. In contrast there are two kinds 
of exchange constants, Jo and J. We assume that 2h 
impurity bonds with Jo randomly distribute among host 
bonds with J in a block of size 26. Then (J^)" 1 is cal- 
culated as [1 + po(J/Jo — 1)]/Js 2 with impurity density 
po = h/b; the randomness does not affect this quantity. 
For (jW) -1 we take an ensemble average since the con- 
tribution of an impurity bond changes whether the impu- 



3 



rity site q is even or odd. We assume that all the possible 
distributions occur in the equal probability and that just 
half of the impurity bonds are on odd sites in the aver- 
age. Hence we have (J^)" 1 = [1 + Po(J/Jq ~ l)]/2 Js. 
The topological angle is given by 9/2-k = 2J^°' / 'J*- 1 ' = s. 
Therefore the bond impurities do not change the gapless 
condition of the homogeneous case p6| . Kawae et al. 
p7[ argued that this model stands in the s = 1 Haldane 
system (CH 3 ) 4 NNi(NO)3 (TMNIN) when nonmagnetic 
impurities Zn 2+ are doped. They observed spin gaps at 
some impurity densities. 

The last example is a site impurity system where impu- 
rity spins are located among host spins. The spin magni- 
tude is so for an impurity spin against s for a host spin. 
Exchange constants are Jo for both sides of an impu- 
rity site and are J otherwise. The number of impurity 
spins in a block is 2h. They are located on the sites of 
{kj \j = 1, 2, • • • , 2h} in the block. Their distribution is 
random in the block and the impurity density is po = h/b. 
In the calculation of s q on Eq. (|ll|), we notice that it be- 
comes s q ^ = X)fc=i(~l) fe+1 s f° r a P ure s P m system. We 
now replace the fcth term in s q by (— l) fc+1 so if an im- 
purity spin locates at the fcth site. Then we have 



■iii 



s(0) 



+ (s-so)£(-l) fej % 



(22) 



where the step function 8(x) is 1 for x > and other- 
wise. Taking the ensemble average with equal weight for 
all possible distributions, we have 



2J (0) i + (i + f)p 



l + 2p 



(23) 



with po = Po{Js/ JqSq — 1). Derivation of this equation 
will be reported elsewhere. When JoSo < J s i the gapless 
condition ( pl|) gives the following results. In the case of 
s = 1, Eq. (|2l|) has no integer solution of I for sq < 2. 
That is, impurities with sq < 2 do not force the gapful 
excitation of a homogeneous s = 1 spin system to be 
gapless p^ , |i"8[ . In the case of s = 1/2, Eq. ( ^H ) has no 
integer solution of / for 1 < so < 5/2. 

In summary, we obtained the NLSM ( |l7| ) for a gen- 
eral antifcrromagnetic Heisenberg spin chain with inho- 
mogeneous spin magnitudes and inhomogeneous nearest- 
neighbor exchange constants arrayed in a finite period. 
We applied this formula to several cases and examined 
the gapless conditions. Since the formula is general, it 
can be applied to various cases which were not treated 
here. Extensions of the present NLSM method to ladder 
chains and two or three dimensional systems are future 
problems. 
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